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Abstract 

The Primitive Equations are a basic model in the study of large scale Oceanic and Atmospheric 
dynamics. These systems form the analytical core of the most advanced General Circulation Models. 
For this reason and due to their challenging nonlinear and anisotropic structure the Primitive Equations 
have recently received considerable attention from the mathematical community. 

On the other hand, in view of the complex multi-scale nature of the earth's climate system, many 
uncertainties appear that should be accounted for in the basic dynamical models of atmospheric and 
oceanic processes. In the climate community stochastic methods have come into extensive use in this 
connection. For this reason there has appeared a need to further develop the foundations of nonlinear 
stochastic partial differential equations in connection with the Primitive Equations and more generally. 

In this work we study a stochastic version of the Primitive Equations. We establish the global ex- 
istence and uniqueness of strong, pathwise solutions for these equations in dimension 3 for the case of 
a nonlinear multiplicative noise. The proof makes use of anisotropic estimates, I/| estimates on the 
pressure and stopping time arguments. 
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1 Introduction 



The Primitive Equations (PEs) are widely considered to be a fundamental model in the study of large scale 
oceanic and atmospheric dynamics. These systems form the analytical core of the most advanced general 
circulation models for the atmosphere (AGCMs) or oceans (OGCMs) or the coupled oceanic-atmospheric 
system (GCMs). Moreover, beyond their considerable significance in applications, the PEs have generated 
much interest from the mathematics community due to their rich nonlinear, nonlocal character and their 
anisotropic structure. 

The physical derivation of the Primitive Equations goes back to the early 20th century [Bje04, Ric07]. It 
is based on a scale analysis that accounts for the relatively constant density of the ocean (in the atmosphere 
the density follows a linear profile) and the contrast between the vertical and horizontal scales (on the order 
of several kilometers vs. thousands of kilometers). We refer the reader to e.g. [Ped82] (or [Ric07]) for further 
physical background for the (deterministic) PEs. 

Hence the Primitive Equations express very fundamental laws of physics and one may wonder what 
the motivations are for introducing uncertainty in "exact" model equations. The introduction of stochastic 
processes in weather and climate prediction is aimed at accounting for a number of uncertainties and errors: 

1. These "exact" models are numerically intractable; they cannot be fully solved with present super 
computers (and will not be for any foreseeable future) . For this reason some sort of statistical averages 
are needed corresponding to the parameterization of the effects of the small unresolved scales. Indeed 
when the full equations are averaged it is common to introduce volumic stochastic forcing terms as a 
partial closure. This is reflected in a growing physics literature on 'stochastic parameterization'. See, 
for example, [Ros77, LL79, MT92, BSLP09, ZFIO]. 

2. The physics of these "exact" models is in fact far from being fully understood. For the atmosphere the 
most dramatic example of physical uncertainty is due to the radiation properties of the air (clouds) 
which produce the source term called Ft in the atmospheric analogue of the equation (1.2d) below, 
that is the energy balance equation. Nowadays these uncertainties due to the radiation properties of 
the clouds are considered to be the most severe source of uncertainty in weather and climate modeling. 
Note that while our presentation focuses on the equations of the ocean (which are mathematically 
slightly simpler than the corresponding atmospheric or coupled oceanic-atmospheric equations) our 
general results and methods in this work apply also to these systems.^ 

Further and related physical background in this connection is the subject of [GHTT]. 

With this backdrop in mind we study in this work a stochastic version of the primitive equations and 
establish the global existence and uniqueness of solutions for a nonlinear, multiplicative white noise. While 
the mathematical theory is now extensive in the deterministic case the stochastic primitive equations are 
only receiving attention very recently. We next review this mathematical background and then conclude the 
introduction by describing our results and outlining the analysis carried out below. 

Previous Mathematical Work 

To the best of our knowledge the mathematical study of the Primitive Equations started in the early 1990s 
with a series of works [LTW92b, LTW92a, LTW93] establishing the existence for all time of weak solutions of 
these equations. Subsequent articles improved these results and obtained existence and uniqueness of more 

^Physical uncertainties on the boundary of the oceanic or coupled oceanic-atmospheric system are also present and are 
responsible for another source of stochasticity in these models. The boundary conditions (1.4) on Vi (which is physically the 
sea surface-air interface) are a simplified version of the more realistic boundary conditions 

av{v - v.) -I-O^v = sv, ariT -n) + d^T = gr- (1.1) 

In the form (1.1) the boundary condition expresses two fundamental aspects of oceanic-atmospheric interaction namely the 
driving force of the wind (the term ) and the heating or cooling of the air by the ocean (the term qt ) . These functions and 
gx are not well known and are estimated by modelers through very rough averages. Thus the uncertainties of these functions 
lead to stochastic PDEs with white noise on the boundary a subject which we intend to pursue in future work. 
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regular (strong) solutions, very similar to the results available for the incompressible Navier-Stokes equations 
(see e.g. [GGMRBOl, HTZ03, PTZ08]). Recently, taking advantage of the fact that the pressure is essentially 
two-dimensional in the PEs (unlike the Navier-Stokes equations) global results for the existence of strong 
solutions of the full three dimensional PEs were established in [CT07] and independently in [Kob06, Kob07]. 
In subsequent work, [KZ07], a different proof was developed which allows one to treat non- rectangular 
domains as well as different, physically realistic, boundary conditions. All of these works make essential use 
of the fact that the pressure terms appearing in the PEs can be shown to be essentially independent of the 
vertical variable z. In the present work we follow an approach closer to that of [KZ07] in that we estimate 
the pressure directly via earlier results for the evolution Stokes equation as in [SvW86]. In any case the 
deterministic mathematical theory for the Primitive equations has now reached an advanced stage and we 
refer the reader to the survey articles [PTZ08, RTT08] for further references and background. 

Notwithstanding these extensive results in the deterministic case, the theory for the stochastic Primitive 
Equations remains underdeveloped. A two dimensional version of the PEs has been studied in a simplified 
form in [EPT07, GHZ08] and more recently in [GHTlla, GHTllb] in the greater generality of physically 
relevant boundary conditions and nonlinear multiplicative noise. While the full three dimensional system 
has been studied in [GII09] following the methods in [CT07], this work covers only the case of additive 
noise. In this case the PEs can be directly studied pathwise via a classical change of variables. In this and a 
companion work [DGHTll] devoted to the local existence of solutions we depart from [GII09] and develop 
different methods which allow us, in particular, to consider a nonlinear multiplicative forcing structure. The 
present work may therefore be seen as the continuation of [DGHTll] which takes us from the local to the 
global existence of solutions. 

In contrast to the primitive equations, the theory of the related stochastic Navier-Stokes equations has 
undergone substantial developments dating back to the 1970's with the initial work [BT73]. In this literature, 
as for the wider literature on stochastic PDEs (see [DPZ92]) two principal notions of solutions have been 
developed. On the one hand Martingale (or probabilistically weak) solutions treat the stochastic elements in 
the problem as an unknown. In this case, one typically establishes compactness in the space of probability 
laws associated to solutions. See e.g. [Vio76, Cru89, CG94. FG95, MR04]. On the other hand Pathwise 
(or probabilistically strong) solutions consider the driving noise as being fixed in advance. See, for example, 
[DPZ96, BFOO, BreOO, BPOO, FR02, DPD03, MR05, Shi07, FR08, GHZ09]. In this and the companion work 
[DGHTll] we concentrate on this latter notion. 

The Governing Equations 

Having described the background and motivations for this article, we now outline its content in more detail. 
The stochastic version of the Primitive Equations that we study below takes the form: 



5tv + (v ■ V)v + w9;,v + — Vp+/k X V - ^JvAv - v^d^^w = + a^{v, T, S)Wi, (1.2a) 

Po 

dzP = -pg, (1.2b) 

V ■w + d^w = {) (1.2c) 

dtT + (v • V)T + wd^T - ^itI^T - vrd^^T = Ft + crriv, T, S)W2, (1.2d) 

dtS + (v • V)S + wd,S - lis AS - vsd,,S = Fs + as(v, T, S)W3, (1.2e) 

p = po{l + Pt{T - Tr) + Ps{S - Sr)). (1.2f) 



Here, U := (v, T, S") = (u,v,T, S), p, p represent the horizontal velocity, temperature, salinity, pressure and 
density of the fluid under consideration; ^v, J^v, Pt, J^t, Ms, vs are (possibly anisotropic) coefficients of the 
eddy and molecular viscosity and the heat and saline diffusivity respectively; / is the Coriolis parameter 
appearing in the antisymmetric term in (1.2a) and accounts for the earth's rotation in the momentum 
equations. The third component of the velocity field, w is a 'diagnostic variable' in that it is determined 
directly from v, the components of the horizontal velocity field (see (2.1b), below). The evolution equations 
(1.2) occurs for (a;i,X2,z) ranging over a cylindrical domain M = A^o x (~^, 0); A^o is an open bounded 
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subset of with smooth boundary dMa. Note that V = {81,82) where 81, 82 are the partial derivatives 
in the (horizontal) xi,X2 directions; A = 9j + 9| is the horizontal Laplace operator. 

The stochastic terms are driven by independent Gaussian white noise processes Wj which are formally 
delta correlated in time. The stochastic terms may be written formally in the expansion 

a^iU)m{t,x) \ / a^{U){x,t)W,^{t) \ 

aTiU)W2it,x) =E 4{U){x,t)W2Ht) , (1.3) 
as{U)W3{t,x) J fe>i V 4iU)ix,t)Wi{t) J 

where the elements are independent white (in time) noise processes. We understand (1.2) in the Ito sense 
but the classical correspondence between the Ito and Stratonovich systems would allow one to treat both 
situations with the analysis herein. We recall the basic mathematical definitions and give precise conditions 
on the operators (Tv, ctt, cts below. 

The boundary dM is partitioned into the top Ti = Aio x {0}, the bottom Ti, = Aio x {—h} and the 
sides Ti = dAio x (— ft,, 0). We denote by uh the outward unit normal to 8J^q. We prescribe the following 
boundary conditions; 

9^v = 0, ^ = 0, 9^T = 0, 9^5* = 0, (1.4) 

on Ti. At the bottom Tb we take 

9^v = 0, w = 0, 8^T = 0, d^S = 0. (1.5) 
Finally for the lateral boundary F/ 

v = 0, d^„T^O, d^^S = 0. (1.6) 

The equations and boundary conditions (1.2), (1.4), (1.5), (1.6) are supplemented by initial conditions for 
V, T and S, that is 

V = vo, T = To, S = So, at t 0. (1.7) 

The Cauchy problem (1.2)-(1.7) given above models regional oceanic flows. We note however that equa- 
tions of a quite similar structure may be given that describe the atmosphere and the coupled oceanic at- 
mospheric system. See e.g. [PTZ08]. The methods developed here could thus be extended to treat these 
systems. 

The manuscript is organized as follows. In the initial Section 2 we set the mathematical background for 
the work defining precisely the notion of pathwise solutions we are studying. Section 3 recalls the results 
in [DGHTll] that guarantee the local existence of solutions. Crucially, these results imply a maximal time 
of existence ^ = C('^); on those samples a; (in the probability space f2) where ^(w) < 00 certain norms 
of the solution U = {v,T,S), in particular the norm, must blow up at £,{uj). Having exhibited these 
preliminaries we next introduce a criteria for global existence based on the uniform control in time oi v G L'^ 
and 8zU G L^. This sets the agenda for the remainder of the paper. In Section 4 we carry out the estimates 
in L^. To this end we introduce a new 'shifted' variable v that satisfies a random PDE that we may analyze 
pathwise. In this way we are able to handle the pressure via the results in [SvW86]. On the other hand a 
number of new terms appear in the nonlinear portion of the equations for v that we must tackle. In Section 5 
we turn to the estimates for 8zU in L^. In this case the pressure disappears when we exhibit the evolution for 
\dzU\'j^2. As such we carry out the estimates in the original variable using stochastic methods: Ito calculus, 
the Burkholder-Davis-Gundy inequality, etc. Finally we include in an appendix further details for various 
technical complements used in the body of the work: the pressure estimates for the Stokes equations after 
[SvW86] and a particular version of the Gronwall lemma that we use to close the estimates in Section 4. 
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2 Mathematical Background and Notational Conventions 



In order to introduce a precise mathematical definition of solutions for the stochastic Primitive Equations 
we begin by rewriting (1.2) in a slightly different form. This formulation will be the basis for all that follows 
below: 



(2.1a) 



1 f° 

dtv + (v • V)v + w{v)d^v + —\/ps - g / (/?tVT + /3s VS") dz 

Po Jz 

+ /k X V - MvAv - iy^d,,w = Fv + av(v,r, S)Wi, 

w{v) = / V • vdz, / V • vdz = 0, (2.1b) 

dtT + (v • V)T + w{yr)d^T - utAT - i^rdzzT ^ Ft + crT(v, T, S)W2, (2.1c) 
dtS + (v • V)5 + w{^^)^,S - /isA5 - vsdzzS = Fs + as(v, T, S)W3. (2.1d) 

Here, as above in (1.2), U := (v, T, S*) = (w, u, T, 5), are the horizontal velocity, temperature and salinity of 
the fluid under consideration, /Zv, J^v, MTj i'T; MSj '^s a-re the coefhcients of the eddy and molecular viscosity 
and the heat and saline diffusivity respectively, / is the Coriolis (rotation) parameter. By integrating 
(1.2b) and making use of the relation (1.2f) we find that the pressure p may be decomposed into a 'surface 
pressure', Ps and some lower order terms that couple the momentum equations to those for the temperature 
and salinity. Crucially, we note that ps does not depend on the vertical variable z. Of course, this system 
(2.1) is supplemented with initial and boundary conditions as given in (1.7) and (1.4)-(1.6) above. 

2.1 Abstract Setting for the Equations 

We next recall the abstract setting for (2.1) (equivalently (1.2)). Note that our presentation and notations 
closely follow the recent survey [PTZ08]. 

Let us first recall some Hilbert spaces associated to (2.1). Define 

H := [{v,T,S) E {L^{M))^ -.V ■ [ vdz ^ in Mo, nn ■ [ vdz^OondMo, f TdM = ( SdM = Q\. 

I J-h J-h JM JM J 

We equip this space with the classical inner product^ which we denote by | • | . Define Ph to be the Leray 
type projection operator from L^{M.Y onto H. For if^(A^)'* we consider the subspace: 

y := ((v,r,S') e (iji(X))^ :V • / vdz = in 7Wo,v = on T/, / TdM = [ SdM = o\. 

I J-h Jm Jm } 

We equip V with the inner product 



(([/,[/«)) 


:=((v,v»))i + 


((T,T«))2 + ((5,5»))3, 


((v,v))i 


:= / (/ivVv 

JM 


■Vv^ +ly^^,^r■^,^r^)dM, 


((r,r«))2 


/ (mtVT 
Jm 


■ Vr« + i^rdzT ■ dzT^) dM, 




Jm 


■VS^ + iysd,S ■d^S*)dM, 



^One sometimes also finds the more general definition (U, C/") := Jj^{v ■ v'd + k^TT^ + KgSS'^)dJV[ with Ky, Kg > fixed 
constants. These parameters kt, i^s a-re useful for the coherence of physical dimensions and for (mathematical) coercivity. Since 
this is not needed here we take kj- = Kg = I. Similar remarks also apply to the space V. 
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and take |1 • |1 = \J ((•, •))■ Note that under these definitions a Poincare type inequality \U\ < c\\U\\ holds 
for aU U €V. We take V(^2) to be the closure of V" n C°°(!M)^ in {H'^{M))''' and equip this space with the 
classical H'^{A4) norm and inner product. 

In the course of the analysis below we shall make estimates involving the individual components of the 
solution U = {u,v,T, S). As such we shall sometime abuse notation and use | ■ | and \\ ■ \\ in the obvious 
way for v = {u,v), T or S. We shall also work with the = LP(A^) norms of these individual components 

of the solution. For p > 1 we denote = {u\u\p-\v\v\p~^) and let |v|lp {Jj^i\u\P + \v\P)dM)^^'' . 
Furthermore, for q,p> 1, we write 




q/p 

{\u\P + \v\P)dz ] dM 



1/9 



Remark 2.1. Assume that v e H^{Ai), with v = on F/ or Jj^ wdAd = 0. As in [PTZ08], an elementary 
calculation that makes use of the Sobolev embedding theorem in M? , reveals that 

Mlu? < c|v|l-lv||■^ (2.2) 

where q >2 and s — I ~ 2/q. This observation will be used on several occasions below. 

The principal linear portion of the equation is defined by'^ 

-/ivAv - v^dzz^^ 
AU = Ph\ -firAT - urdzzT I , for any U = (v, T, S) £ D{A) 



where: 



-Hs^S - vsdzzS J 



D{A) = {U^ (v, T) e V(2) :a,v = = dzS = on F,, 

d^„T = d^„S - on Vi,dzV ^ dzT - 9,5 = on FJ. 



(2.3) 



We observe that A is self adjoint, positive definite. Indeed, by integration by parts and using the boundary 
conditions imposed by (2.3), we see that {AU,U^) {{U,0)) for aU U,0 G D{A) and thus, by density, 
for all U, G V. Note also that, due to regularity results for the Stokes problem of Geophysical Fluid 
Dynamics, \AU\ ^ \U\h2. See [Zia95] and also [PTZ08]. 

We next turn to the quadratically nonlinear terms appearing in (2.1). Noting that there is no momentum 
equation for w in (2.1) and in accordance with (2.1b) we define the diagnostic function; 

w(f/) = w(v) /" V-vdz, U = {v,T,S)&V. (2.4) 

J Z 

Take, for U,0 e D{A): 



(v • V)v«^ 

Bi(C/, C/«) Ph\ (v • V)r« I , B2(C/, C/«) := Ph 
,(v-V)5», 



/ w(v)9,v«^ 

wiw)dzTi I . (2.5) 
\wMdzS^, 



We let B{U, U^) := Si(C/, U^) + B2{U, U^), and often write B{U) = B{U, U). As in [PTZ08] one may show 
that B is well defined as an element in H for any U,U^ G or V(2). In addition to the properties of 

B appearing in [PTZ08] we have the following additional bounds, which are established with anisotropic 
estimates along the same lines (cf. Remark 2.1). 



^In comparison to previous works, such as [PTZ08], we do not include all of the terms due to the pressure in the definition 
of A. Such elements destroy the symmetry of A and arc therefore relegated to a lower order term Ap. See (2.7). 



6 



Lemma 2.1. Suppose that U,U^ £ D{A) and that e H. Then 



\{BiU,U^),U')\ <c(|vU4|lC/«||l/4|^C/»|3/4|c;b| ^ ||^||l/2|^|l/2|^^^S|l/2||g^^S||l/2|f^b|)^ (2.6) 



1(2) 

For the second component of the pressure in (2.1a) we take 

-gJ°{l3TVT + (3sVS)dz 



ApU^Ph\ \, U&V. (2.7) 



We capture the Coriohs (rotational) forcing according to 



Finally we set 



/k X V 

EU = Ph \ ' \ , U eH. (2.8) 



(2.9) 




We shall assume throughout this work that: 

FeL^in;Ll,i[0,^);L\M))). (2.10) 

We finally give a precise definition for the stochastic terms appearing in (2.1) (i.e. (1.2)). For this 
purpose let us briefly recall some aspects of the theory of the infinite dimensional Ito integration. As we 
are studying pathwise solutions of (1.2) (see Definitions 2.1, 2.2 below) we shall fix throughout this work a 
single stochastic basis S := {il,J^,¥, {J-t}t>o, W). Here is a cylindrical brownian motion defined on an 
auxiliary Hilbert space il and adapted to the filtration {J-t}t>o- By picking a complete orthonormal basis 
{efc}fe>i for if, W may be written as the formal sum W{t,uj) = J2k>i SkWk{t,oj) where the elements Wk are 
an independent sequence of ID standard Brownian motions. 

Consider another separable Hilbert space X and let L2{il,X) = {R £ C{!d,X) : X^fcl-^'^feP ^ 
that is the collection of Hilbert-Schmidt operators from il into X. Given an X-valued predictable^ process 
G G L^{il;Lf^^ {[0, go), L2{il, X))) one may define the (Ito) stochastic integral 



Mt := / GdW = V / GkdWk, 
Jo J. Jo 



as an element in Aix, that is the space of all X- valued square integrable martingales (see [DPZ92] or [PR07]); 
here Gk = Gek- The process {Mt}t>o has many desirable properties. Most notably for the analysis here, 
the Burkholder-Davis-Gundy inequality holds which in the present context takes the form. 



E sup 

\ t6[0,T] 



GdW 




\ r/2 

\G\l,in.x)dt] , (2.11) 



valid for any r > 1. Here c is an absolute constant depending only on r. 



*Let <I> = f2 X [0, oo) and take Q to be the cr-algebra generated by sets of the form 

{s,t]x F, < s <t < oo,F e Ts; {0} x F, F &To. 

Recall that an X valued process U is called predictable (with respect to the stochastic basis cS) if it is measurable from {^,Q) 
into {X, B{X)), B(X) being the family of Borel sets of X. 
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Given any pair of Banach spaces X and y we denote by Bndu{X,y), the collection of all continuous 
mappings "if : [0,oo) x X ^ y such that 

\\-9{x,t)\\y <c{l+\\x\\x), xeX,t>0, 

where the numerical constant c may be chosen independently of t. If, in addition, 

- 'if{y,t)\\y < c\\x ~ y\\x, x,yeX,t>0, 

we say is in Lipu{X, y). We define 

/ ^v(C/) \ 

a{{v,T,S)) = a{U)^PH\ MU) \, U e H, (2.12) 

V '^siU) J 

and assume that a : [0, oo) x iJ — > L2{il, H) with 

a G L^pu{H,L2{ii,H))f^Lipu{V,L2{'d,V))f^Bndu{V,L2{'d,D{A))). (2.13) 

Under the assumption (2.13) on cr, the stochastic integral t ^ a{U)dW may be shown to be well defined 
(in the Ito sense), taking values in H whenever U S L'^{fl,Lfg^{[0,oo);H)) and is predictable. 

Remark 2.2. The condition (2.13) may be shown to cover a wide class of examples, including but not limited 
to the classical cases of additive and linear multiplicative noise, projections of the solution in any direction, 
and directional forcings of Lipschitz functionals of the solution. See [GHTT] for further details and physical 
motivations. 

We note that the final condition a G -Bnc?u(V, i2(ii, -D(j4))) in fact may be weakened slightly to a £ 
Bndu{D{A), L2{ii, D{A))) for the proof of local existence. However for the global existence we need this 
stronger condition in order to carry out the estimates which appear in Section 4- In fact the stronger 
condition used in this section is an artifact of a 'pathwise ' approach which further complicates the nonlinear 
term. This pathwise approach is necessitated by the appearance of certain terms explicitly involving the 
pressure Ps . 

Collecting the operators defined above we reformulate (2.1) (equivalently, (1.2)) as the following abstract 
evolution system 

dU + {AU + B{U)+ApU + EU)dt = Fdt + a{U)dW, U{0) = Uo. (2.14) 

We recall the following pathwise (or probabilistically strong) notions of local and global existence for this 
system. 

Definition 2.1 (Local Solutions). Let S = (51, J-, {J^t}t>o, W) be a fixed stochastic basis and suppose that 
Uq is a V valued, J-q measurable random variable with E||[/o||^ < oo. Assume that F satisfies (2.10) and 
that (2.13) holds for a . 

(i) A pair (JJ,t) is a a local strong pathwise solution of (2.14) (i.e. of (1.2)j if t is a strictly positive 
stopping time and U{- At) is an Tt-adapted process in V so that (relative to the fixed basis S) 

[/(•At) GL2(fl;C([0,oo);y)), f/]lt<, G ^^(r!; LL([0, oo); 7^(A))); (2.15) 

and, for every t > 0, 

U{tAT)+ {AU + B{U) + ApU + EU)dt' = Uo + Fdt' + a{U)dW, (2.16) 
Jo Jo Jo 

with equality understood in H . 



8 



(a) Strong pathwise solutions of (2.14) are said to he (pathwise) unique up to a stopping time r > if given 
any pair of solutions {U^,t), (U'^,t) which coincide at t = on the event Q = {[/^(O) = C/^(0)} C il^ 
then 

P(ljj(C/Ht At) - U^itAr)) = 0;Vi > O) = 1. 

Remark 2.3. For a local solution {U, t) consider the increasing function 

0(cr) := sup \\U{t)f {defined for a e[0,T]). 
te[o,(T] 

We infer from (2.15) that 4>{t) < oo, almost surely. Of course the question arrises as to whether the solution 
U and therefore (j) can be extended beyond r. We will call ^ the supremum of the t's such that {U,t) is a 
local solution of (2.14). A priori, ^ might be finite or infinite and we are lead to the following notions of 
maximal and global solutions. 

Definition 2.2 (Maximal and Global Solutions). 

(i) Suppose that {t„}„>i is a an increasing sequence of stopping times converging to a (possibly infinite) 
stopping time ^ and assume that U is a predictable process in H. We say that the pair (U,^) is a 
maximal pathwise strong solution if (U,Tn) is a local strong pathwise solution for each n and 

</,(^) = sup \\Uf = oo (2.17) 
te[o,C] 

almost surely on the set < oo}. 

(a) If (U,^) is a maximal pathwise strong solution and ^ = oo a.s. then we say that the solution is global. 
Note that (j){oo) can be finite or infinite in this case. 

With these foundations in place we finally state, in precise terms, the main result of the work, the global 
well-posedness of (1.2) in the class of pathwise, strong solutions: 

Theorem 2.1. Suppose that Uq G L'^{rt,V) and is J-q measurable. Assume that F satisfies (2.10) and that 
(2.13) holds for a. Then, there exists a unique global pathwise solution U of (2.14) with U(0) — Uq. 

The proof of this Theorem is given below at the end of Section 3. It follows directly from the combination 
of Theorems 3.1, 3.2 and Propositions 4.1, 5.1. 



3 Local Existence and A Criteria For Global Existence 

As with the deterministic case, the first step in the analysis is to establish the local existence of solutions 
all the way to a 'maximal' existence time .J > 0. Recall that if £^{uj) < oo for w S fi, then the norm of 
the solution must blow up at this maximal time as expressed by (2.17). With this in hand, we next show in 
Theorem 3.2 that the strong norms of solutions may be controlled by a combination of the L'^ norm of the 
momentum v and by the vertical gradient of the entire solution U. The results in Theorem 3.2 establish a 
sufficient condition for Theorem 2.1 to be valid and thus set the program for the rest of the article. 

The first step, establishing the existence of a maximal pathwise solution, is given by Theorem 3.1 imme- 
diately below. The proof of this result follows immediately from abstract results in [DGHTll, Theorem 2.1 
and Section 6]. 

Tiieorem 3.1. Suppose that (2.10) and (2.13) hold for F and a respectively. Given any J-q -measurable 
Uq G L'^{Q,,V), there exists a unique maximal pathwise solution (U,^) of (2.14) with U{0) = Uq. If, for 
some p > 2, [/(O) G LP{Q, V) then, for any t > 0, we have the bounds 

e( sup \U\P+ [ ^ \\U\\^\U\P-^dt'] < oo, (3.1) 
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eU \\U\\^dt'\ <oo. (3.2) 

Remark 3.1. The local existence of solutions for a (nonlinear) multiplicative noise is more challenging than 
the additive case. In this latter situation where the noise terms have no state dependence, a simple change of 
variables which 'subtracts off ' the noise can be used to transform the governing equations into a random PDE 
which can be treated using the methods of ordinary calculus. Here, crucially, the probabilistic dependence, 
given as oj € may be essentially treated as a parameter in the problem. 

In the general case of nonlinear multiplicative noise no such transformation is available. As a consequence 
stochastic integral terms must be estimated directly and as such the classical compactness methods break down 
since the problem may no longer be handled with the 'pathwise' approach typically used in the additive case. 
More involved probabilistic machinery is therefore needed. The proof of Theorem 3.1 appearing in [DGHTll] 
is based on an elementary but powerful characterization of convergence in probability given in [GK96]. With 
this characterization we were able to establish a Yamada- Watanabe type result namely: 'pathwise ' solutions 
follow from the existence of Martingale solutions when we have pathwise uniqueness for solutions ( see Defi- 
nition 11, (ii)). A different method based on the Cauchy convergence of the Galerkin solutions associated to 
the basic SPDE is developed in [GIIZ09J for the Navier-Stokes equations in dimensions D = 2,3 and for the 
2D Primitive equations in [GHTlla, GHTllb]. 

For any give A' > we define the following stopping times which are used here and below: 
r^:= inf ( sup \U\' + r\\\Ur + \F\l,)dt' > k] 

= supi sup \U\^+ [ \\\Uf + \F\l,)dt' < k\ , (3.3) 
t>o lt'e[o,tA{] Jo I 



and also 



r]^^ := inf < sup |v|^4 > K > = sup < sup |v|^4 < K 



t'G[o,tA5] I t>o lt'e[o,tA{] 



4'^:=inf| sup |9.C/|'+ / ^la.C/f >i^l =sup I sup \d,Uf + [ 
t>o |^t'G[o,tA4] JO J t>Q [t'e[o,tA{] Jo 



tAC ^ 



rK:^4^ /\r^K^. (3.4) 

Note here that we employ notation typically used in optimization (see e.g. [ET99, Roc70, Pol71, Pol97]) 
and we follow the convention that the infimum of an empty set is +oo. By Theorem 3.1 and the standing 
assumption (2.10), it is direct to infer that limx-|-oo = oo. We establish below, by gathering together the 
conclusions in Propositions 4.1, 5.1, that Iuhk^ooTk = oo. This, as the next result shows, is a sufficient 
condition for the global existence of solutions of (1.2). 

We observe and emphasize that the definition of tk does not imply that tk < S, a.s. Indeed, it could be 
possible that Ti^ = oo if the sets defining tk are empty. Rather, that tk < will be proven as part of the 
next theorem. 

Theorem 3.2. The hypotheses are the same as in Theorem 3.1 and we consider the maximal strong pathwise 
solution {U,£,) = ((v,r, S*),^) of (1.2) define in this theorem. For any deterministic constant K, consider 
the stopping times tk defined according to (3.4). Then, for any deterministic time t > Q and any K > Q, 

/ pt/\TKA£, \ 

E sup \\Uf+ \AUfdt' ]<oo. (3.5) 

\t'<£[0,tATK/\i] JO / 
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Furthermore tk < ^ almost surely and consequently if we show (by other means) that limA'^ ta' = +00 then 
^ = cx) a.s., that is they solution U is global in the sense of Definition 2.2, (ii). 

Proof. With an application of the Ito formula we deduce the following evolution equation for ||C/|[^: 

d\\Uf + 2\AU\^dt = -2{B{U) + ApU + EU + F, AU)dt + \\a{U)\\l^(^^ y-^dt + 2{A^^^a{U), A^/^U)dW. 

Fix any < Ta < Tb < A ^ A i. Integrating in time, taking a suprcmum over the interval [tq, ti,] and then 
the expected value of the resulting expression yields the estimate 



E( sup \\U 

\Ta<t'<Tb 



2 j " \AU\^dt' 

<E||C7(Ta)f +cE f ' \{BiU),AU)\dt + cE I {\{ApU + EU + F, AU)\ + \\a{U)\\l^^^ y))dt 



cm sup 

Ta<t'<Tt 



{A^/^a{U),A^/^U)dW 



■.=E\\U{Ta) 



J1 + J2 + h 



We address the terms Jk in reverse order. For J3 using the Burkholder-Davis-Gundy inequality, (2.11), and 
the assumption in (2.13) that a G Lipu{V, L2{'d, V)) we infer 



1/2 



J3<cE\^j {A^/^a{U),A^/^U}^dtj < cffi 
< ^E sup ||C/f + cE / {l + \\Uf)dt. 

2 t'elTa.Tt] JTa 



HU)\\ 



L2iii,V) 



lUW^dt 



1/2 



For J2 we use again that a e LipuiY, L2{iX, ^)) '^^'^ make direct estimates for Ap and E; see respectively, 
(2.7), (2.8) above. We infer 

\{ApU + EU + F,AU)\ + |k(C/)|li^(u,^) < ^lAUf + c(l + \\Uf + \Ff). 
On the other hand, (2.6) implies that 

\{BiU),AU)\ < \w\l4U\\'/MU\'/' + \\U\\'/'\dM^'\\d^U\\'/'\AU\'/') 
<l\AUf+ciW\l. + mf\mf)\\Uf. 
This takes care of Ji. Combining these three estimates, we conclude that 



E( sup \\U\\ 

\Ta<t'<Ti, 



\AU\^dt' 

<cE[\\u{Ta)W'+ r\Ml. + \d,u\-'\\d,ur + i)\\urdt+ I \F\-'dt 



(3.6) 



We now apply a stochastic version of the Gronwall lemma to (3.5) from [GHZ09, Lemma 5.3] which is 
recalled below in the Appendix as Proposition 6.3. Note carefully that the constants appearing in the above 
equation do not depend on Ta,rb. By taking X := \\U\\^, Y := \AU\^, Z := \F\^, R := \v\l^ + \dM\^\\d^U\\^ + l 
and finally t := t A tk ^ ^ in Proposition 6.3 we therefore derive (3.5) from (3.6). 

We now argue by contradiction to show that tk < ^. If. on a set of non-zero measure, tk > S, then using 
that {tk > ^} = U(>o{Tif A i > ^} we would infer the existence of a deterministic time t > such that 
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{tk At > £,} is oi non-zero measure. By the definition of ^ (cf. (2.17)), we would therefore infer that for this 
t > 0, 

sup ||?7f + / \AU\''dt'> sup \\U\\^ = oo 

t'e[o,tATK-AC] "'0 t'e[o,C] 

over {tk At > £,}. Since {tk A i > ^} is a set of non-zero measure we have a contradiction witli (3.5). With 
this, the proof of Theorem 3.2 is therefore complete. □ 

We conclude this section by explaining in precise terms how we will use Theorem 3.2 to prove Theorem 2.1: 

Proof of Theorem 2.1. Define the stopping times rj^', t^^\ and take tk Arj^' as in (3.4). According 

to Theorem 3.2. the desired result will follow once it is shown that lim/f-|-oo tk ~ oo a.s. In Proposition 4.1 
it is shown that limif-|-oo t^"^ — oo and from Proposition 5.1 we infer that limi^ioo = oo. The proof is 
therefore complete. □ 



4 Estimates in L^(A^) 

In this section we establish that Iuhk-^oo t'k'' = +oo where rj^'' is defined as in (3.4) above. This amounts to 
showing that the horizontal component of the velocity field, v, of the solution U = (v, T, S) is almost surely 
finite in L'^L^ ^, at least up to the maximal time of existence, ^. 
More precisely we will prove the following: 

Proposition 4.1. Suppose that the conditions given in Theorem 2.1 hold and consider the resulting maximal 
strong solution {U,£,) = ((v,T, 5*),^) of (1.2). Then, for any t > 0, 

sup |v(t')|i4 < oo a.s. (4.1) 

t'e[o,tA5] 

Moreover, defining, for K > the stopping time r|3'' as in (3.4) we have, up to a set of measure zero, that 
limiffoo t'^k = oo- 

Proposition 4.1 will be an easy consequence of Propositions 4.2, 4.3 below. The proof relies on a decom- 
position of the solution U of (2.16) into a sum of two components: JJ = (v, T, S) and tj = (v, T,S); U solves 
a linear stochastic system, (4.2), whereas U is the solution of the random PDE, (4.4). The later system 
may be studied 'pathwise', that is pointwise in fl. We are therefore able to treat certain terms involving the 
pressure explicitly by means of Proposition 6.1 (cf. [KZ07], [SvW86]). In any case, we arc able to establish 
control on the norms ^ of v and v as in (4.1) in Propositions 4.2, 4.3 from which Proposition 4.1 then 
follows immediately. 



The Decomposition 

We describe next the decomposition and recall some basic properties of the linear portion of this splitting. 
Consider the following linear system, an infinite dimensional Ornstein-Uhlenbeck process: 

dU + AUdt ^ lt<iCT(U)dW, U{0) = 0. (4.2) 

We underline here that the element {U, ^) appearing in the right hand side is the unique maximal solution 
of (1.2) corresponding to the given initial condition Uq. 

Proposition 4.2. There exists a unique global, pathwise strong solution U = (v,r, S*) of (4.2) with U 
e L'^{il;C{[0,oo); D{A))). In particular, for any t>0, 

sup |v(t')|i4 < oo a.s. (4.3) 

t'e[o,tAC] 
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This result is proven exactly as in [GHTllb, Lemma 4.1]. See also [DPZ92] for further generalities 
concerning infinite dimensional linear stochastic systems of the form (4.2). 

Remark 4.1. The final condition in (2.13), a G Bnd{V, L2{iX, D{A))) is imposed so that U is guaranteed to 
evolve continuously in D{A). This stronger conditions is needed in several estimates below (see e.g. (4.11), 

(4.15);. 

We next define U — U — U. By subtracting (4.2) from (2.14) wc find that U satisfies 



d 

U{0) = C/(0) = Uo. 



— U + AU + B{U + U)+E{U + + Ap{U + U) = F, 



The equation for the first component v of the solution U ~ (v, T, S) is given by 

dtV + ((V + V) • V)(V + V) + W(V + V)9^(V + V) + —VPs + /k X (V + V) - ^vAv - I^v^^zV 

Po 



+ g ^ (/3t V(f + f)+ f3sV{S + 5)) dz. 



(4.5) 



Note that the diagnostic function, w{-) is defined as above in (2.4). As, with v, we have /^^ V • irdz = 0. Of 
course, (4.5) is supplemented with boundary conditions as in (1.4)-(1.6). We will hereafter prove: 

Proposition 4.3. Define (v,T, S) = U :=U — U . Then for any, for any t > 0, v satisfies 

sup |v(t')|i4 < oo a.s. (4-6) 

t'6[0,tA5] 

Proof. Note that (4.5) does not have any white noise terms. Multiplying (4.5) by v'^ = {^1,^2) and then 
integrating over the domain M. leads to the following system describing the time evolution of |v|^4: 

J A: "'■^ A: ^ 

= - V / {V]+ Vj)djVkvldM - V / w(v + w)d^VkvldM 

-Y—f djPsvUM + q I I (/3TV{f + f) + /3sV{S + S))dz-v'^dM 
, Po Jm JmJz ^ ^ 



(4.7) 



- / / k X (v + v) . <f^dM + I Fv • v^dTW 



:=Ji + J2 + J3 + J4 + J5 + Je- 
Note that we have made use of the cancelation property: 

(v« • Vvv'' + w(v»)a,v • v'')d7W = 0, (4.8) 

M 

which holds for any ?■ > 1 and any vector fields over A4, v, that arc sufficiently smooth and satisfy the 
boundary conditions (1.4), (1.5), (1.6). Notice moreover that for some constant k > (depending only on 
jUv7 i^v) we have that 

^\V,irr<^Y.f id.vuY^ldM+'^Y.f id.^^f^ldM. (4-9) 
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where V3 = {di,d2,dz) is the full three dimensional gradient operator. 



Estimates for the Nonlinear Terms Ji, J2. 

We first provide estimates for the inertial terms Ji, J2 in (4.7). For Ji we estimate: 

k,j •'M Jm \Jm J 

^^l^t^lE (/^ i^'f^'dM) < c\AU\ (|v^|^/^|V3(v^)|i/^) 

<c|^J7||v|i.|V3(v2)|V2<c|AL/r/3|v|i.+..|V3(v2)p. 

Here we have used the Sobolev embeddings in M'^ of into and of H^^'^ into 

For Ji we also need 

to estimate 



^ / VjdjVkvldM 



< c|v|i^|Vv|i4|v3|^,/3 < c\AUWl + 



In summary, 



|Ji|<c(l + |^J7n(l + |v|10 + «|V3(v2)|^ 



The estimates for J2 make use of a preliminary integration by parts (cf. (2.1b)) 



= J2 [ I {vj + i'j) dz d,,VkvldM + [ [ 

j f. J M J z J 1^ J M J z 



{vj + Vj) dz dzVkdjVkVkdA4 



■ — J2,l + J2.2 + J2,3 + J2,4- 



(4.10) 



For J2,i we use the embedding of i/"^/* into L'* in and find, 



\J2A\<c\v\Ls\Vdz^\\v'\LS,s<c\AU\\^\ls<c\AU\Y(^J^ 



i^) ) 



^\ (l/4)-2 



< 



c\AU\ (|v^|V4|V3(v^)|^/^) < c\AU\'\m. + f |V3(v2)p. 
For J2,2, by making use of Agmon's inequality and that H^^^ is embedded into we infer: 



\J2,2\ <|vUoo|V5.v||v3| < c|v|ioo|V9,v||v2|^/2 < c\AU\'' (\^'\'^'\y3i^'')\'^'Y^' < c\AU\'\v\%^\V^i^r')\^^' 



<c\AU\''/'\v\lT + f |V3(v^)|^ < c(l + \AU\'){1 + |v|i.) + f |V3(v^)r. 
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J2,3 and J2,4 seem to require 'anisotropic' type estimates: 



J2,3 + J2,4\ / {Vj+Vj)dz 



\dzVk\L<i\djVkVk\L^, \vk\L^^dMo 



< / \vj + Vjl^ldzVklLildji'kVklLilvklL'idMo 

Note that we made use of Remark 2.1, (2.2) with = 12 in order to estimate j-Oj+Wj [^2^12 = |tij|^2^i2 < c||J7||. 
Summarizing the above estimates and taking advantage of the observation (4.9) we conclude that 

IJ2I <c(l + |lC/f)(l + |^[/n(l + |v|iO+Aiv^ / {d,VkfvldM + v^Y. I {dzVufvldM. un) 



Estimates for the Pressure Term J3 

We next attend to the term J3 . Using the fact that the pressure term, ps , depends only on the horizontal 
variable x we find 



IJ3 



V — / djPs [ v^dzdMo 
, PO J Mo J-h ^ 





.0 


< c|Vp,, 1^4/3 




J-h 



Li 



By applying the Sobolev embedding in Z? = 2 that Vt^^/'^'i is embedded in L^{M.o) we have 



w dz 



-h 



j.k 
j.k 



dj I vidz 

-h 



^^E 



Ma \J-h 



{djVkYv'idz 



{djVkVk)vkdz 



-h 



r 4/3 



\ 2/3 / ^0 

a- 2 



2/3 

vldz \ dMr 



-h 



3/4 



<c|v|l4 Y.(( {d.VkfvldM 



1/2 



Combining these estimates and using Young's inequality we conclude that: 

IJ3I < c|Vp.|2./3|v|i. +/iv^ / (djVkfvldM. 



(4.12) 



Remark 4.2. Note that it is at precisely this stage in the estimates that we see the crucial role played by 
the two dimensional spatial dependence of the surface pressure Ps ■ This insight concerning the importance of 
the lower dimensional pressure is the key to the recent breakthroughs for global existence in [CT07, Kob06, 
Kob07, KZ07]. In the present study we follow the approach in [KZ07] which treats the pressure explicitly via 
bounds on the pressure for the Stokes equation found in [SvW86] and reproduced below in Proposition 6.1 for 
the convenience of the reader. 
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Ill order to estimate the term |Vps|^4/3 appearing in (4.12) we now apply Proposition 6.1 below. To this 
end we average (4.5) in the vertical direction. Define the operator 2l(v) = i J^^^vdz and let q = (Oo2l(v). 
By applying 21 to (4.5) we find that q satisfies the following Stokes system over Aio C M^: 

9tq - AivAq + Vp, = G(C/), 
V • q = 0, q\0Ma = 0: 



where, for v = v + v, 

GiU) = G(v, T, S) = Gi{U) + G2{U) 



-Pq% ((v • V)v + (V ■ v)v) + poSl {^9 {Pt^T + PsVS) dz - fk x v + 



(4.13) 



Note that, using (1.2c) and taking into account the boundary conditions (1.4), (1.5), 2t(w(v)9zv) = 2t((V ■ 
v)v) and 2l(9;,^v) = 0. 

Making use of Proposition 6.1 below, we find, for any < t < oo, and for any pair of stopping times Ta, 
Tb, with 0<Ta<Tii<tA^ that 

NPsll.J.dt' < C (^MTaW+ £ \G{U)\l,jsdt' 

Observe that l|q(Ta)|[ < c||v(Ta)|| and that 



\G2{U)\l,,sdt' < c / \G2{U)\^dt' < c I (||[/f + \F\^)dt'. 



On the other hand. 



^ ' \Gi{U)\l,jsdt' <cj^ " |(v . V)v + (V ■ ^)v\l,,,dt' 

(j |v|^/3|Vv|^/3dA^j dt'<cj \^\l4vf dt'. 



In conclusion: 



£ mil^jsdt' < c (^||v(rjf + £{{1 + Mh)\\Uf + \F\')dt'^ . (4.14) 
Finally, combining (4.14) and (4.12) we find, for any pair of stopping times < Tq < Tt < t A ^, 



'\J3\dt'<l sup |v|i4 + Mv / ' V / {d.VkfvldMdf 



\ (4-15) 



+ c l|v(rjf + (l+ sup (|v|i. + |v|i.) / {\\Ur + \F\^)df 

\ t'e[Ta,T6] Jt^ J 

Estimates for the Lower Order Terms J4, J5, Jq 

To estimate the term J4 we use the embedding of into H^^^ in 2D and we observe that 

/ /• \ 1/3-3/2 rt/n 

lh\ < c\\U\\ J2 [JjilfdMj < c\\U\\ (|v2|V2|V3V^|V2) 

<«|V3V^|^ + c(l + ||[/|r)(l + |v|i.). 
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The estimate for the Coriohs term J5 is direct 

|J5|<c|vU.|v|i. <c|lC/||(l+|v|i.). (4.17) 

Finally for Jg wc find 

iJel < |i^v|L*|v|i. < c{\F\l, + l)(|v|i, + 1). (4.18) 

Integrating (4.7) in time from Ta to Tb and then putting together the estimates (4.9) (4.10), (4.11), (4.15), 
(4.16), (4.17), (4.18) we find that for any pair of stopping times t^, with 0<Ta<Ti,<tA£, 

sup |v|i4 <c||v(T„)f + c(l+ sup |v|i4+ sup \ir\iA f ' H{t')dt', 

t'6[ra,T6] y t'£[Ta,Tb] t'e[0,tA?] J J 

where 

H{t') := (1 + ||C/(t')in(l + \AU{t')f) + \F{t')\l,. 

We finally apply the version of Gronwall's lemma in Proposition 6.2 below with p = 2, X{t') = |v(t')|^4, 
f{t') = c||v(t')|p, g{t') = cH{t') and h{t') = c(l + sup^/^jQ j^^j \\r\'^i)H{t'). Observe that, as a consequence 
of Theorem 3.1 and Proposition 4.2 it is clear that U ~U — tj <E C{[0,t A^); V). It therefore follows directly 
that X e C{[0,t A ^)) and that / is continuous at t = 0. Also, due to Theorem 3.1, Proposition 4.2 and the 
standing assumption (2.10), we easily infer that f,g,h G L^{[0,t A^]). Thus by. Proposition 6.2, we conclude 
that (4.6) holds and Proposition 4.3 is now proved. □ 

Proof of Proposition 4-1 (conclusion). Combining (4.6) with (4.3) wc infer (4.1). By the definition of rj^'', 
(3.4), it follows directly from (4.1) that lim/^-f-oo t]^'' = 00, this completes the proof of Proposition 4.1. □ 

5 Vertical Gradient Estimates 

In this section we will show that limx-foo = +00, with rj^ defined in (3.4). Recall that, in view of 
proving Theorem 3.2 and hence establishing the main result of the article, we need to show that liniK'foo tk ~ 
+00, where tk — t^^^ A rj^^ is given by (3.4). Since we have already shown in the previous section that 
limiffoo = +00 we now need to show that limx^ which essentially amounts to deriving 

suitable bounds for dzU in L^{Ai). 

Proposition 5.1. Assume the conditions set out in Theorem 2.1 and suppose that (U,^) is the corresponding 
strong pathwise solution of (1.2). Then, for any t > 0, 

sup + / WdMfdt' < GO. (5.1) 

t'sKAt] Jo 

Hence, defining rj^'' according to (3.4), limA'-foo t]^' = 00. 

The boimd (5.1) and hence Proposition 5.1 follows directly from (5.4) and (5.18) which are established 
in Corollary 5.1 and Proposition 5.3 respectively below. For the analysis we make use of several additional 
stopping times: 

sup \dMtT+ \\dMt')f >K} ^snp I sup |a.v(t')|' + / ||a.v(t')f < A' 
6[o,tAC] Jo J t>o [t'e[o,tA5] Jo J 

= inf<; sup (|r(t')li4 + |5(t')li4)>if| =sup| sup {\T{t')\l^ + \S{t')\U<K], 
t'e[o,tA4] J t>o [t'e[o,tA?] J 

Ar^^^ArW. (5.2) 





inf I 
t>0 1 


T,S 

Tk 


:= inf I 
t>o 1 




— T^'^ 
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Here is defined as in (3.4). We first tackle dz^ in L^. From these estimates along with those carried 
out above for v in L'^{Ai) wc infer that \im.K-ioo = oo. We next would like to estimate d^T and dzS 
in LP'^M.). However, in order to be able to treat the nonlinear terms which arise (see, for example, (5.20) 
below) we first must estimate T and S in L^(7W). This is the reason for the second collection of stopping 
times, , are defined in (5.2). In this case, in contrast to the momentum equations above, there are no 
pressure terms which means that we can perform these estimates in the original variables by means of 
the Ito lemma. With these estimates we deduce that Vlthk^oo ctk = oo which allows us to finally establish 
the desired bounds for dzT and dzS, (5.18), in Proposition 5.3. Having now established (5.1), we thus 

(2) 

immediately infer that limA'^oo ''"^ ^ desired for Theorem 3.2. 



5.1 Estimates for the Momentum Equations 

Proposition 5.2. Assume the conditions set out in Theorem 2.1 and suppose that (C/, ^) is the corresponding 
strong pathwise solution of (1.2). Define r^'' as in (3.4). Then, for any t > 0, 

( /-rlc^ACAt \ 

E sup |9,vp + / Wdz^fdt' < oo. (5.3) 

yt'G[0,rj^''ACAt] "'0 J 

From this proposition, proven immediately below, wc obtain the following 

Corollary 5.1. As.sume that the conditions set out in Theorem 2.1 hold and that (U,^) is the resulting 
strong, pathwise solution of (1.2). Then for, any t > 0, 



sup la^vp + / \\dzv\\^dt' < oo a.s. 
t'effAtl Jq 



'eKAt] 

Moreover, defining the stopping times t^'^ according to (5.2) we have that lim;f-|-oo ~ 

Proof of Proposition 5.2: We apply dz to (2.1a) to infer an evolution equation for dzV. 

ddzv + {dz[{v ■ V)v + w{v)^z^r] + gdzW^T + /3s V5] + /k x (a,v)) dt 
- in^AdzV + v^dzzz^) dt = dzF^dt + a^crv(v, T, S)dW. 



(5.4) 



(5.5) 



Note that, since Ps is independent of z, the pressure term does not appear above. After an application of 
the Ito formula to this system we find 

d|a,v|2 + 2||a,vf = - 2 / (a,[(v • V)v + i(;(v)a,v] +59,[/3tVT + /35V5]) • dz^dMdt 

+ 2{dzF^, dz^)dt + |9,av(v, r, 5)|i^(,,^)di + 2(9,av(v, T, 5), dz^)dW ^^'^^ 
= ( Jf + J| + J| + J| + J|)dt + J^dW. 

Here we have made use of the boundary conditions (1.4), (1.5), (1.6) in order to infer (— (/ivAS^v + 
j/ycJzzzv), 9zv) = ||(9zVp. Note also that there is a cancelation in the 'Coriolis term' [k x dz'v,dz'v) ~ 0. 
We begin by treating the nonlinear terms. For J^, an integration by parts reveals 



Jf = - / d^VjdjVkdzVkdM - / vjdzdjVkdzVkdM 

= djdzVjVkd^VkdM + dzVjVkdjdzVkdM I v^d^^djVkdzVkdM. 
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Hence with direct estimates using Holder's inequality and the Sobolev embedding of H^/^ into L^{M.) 

\Jl\ <c||5,v|||vU4|9,vU4 < c||5,v||V4|v|^,|a,v|i/4 < c|v|i.|9,vp + \\\dM? ■ (5-7) 
For the second portion, J|, of the nonlinear term we have: 

J| = / w{w)d,w ■ d,,wdM = ]- I w{w)d, {d,y ■ d,v) dM 

2 Jm „■ Jm 



and hence, as for we may estimate 



kll<c|v|i4|a,vp + -||a.v|| 



Direct estimates and using the assumption (2.13) lead to 



I J|| + 1 J|| + 1 J5I < c(i + \\ur + \F\') + -iia.vf . 



(5.8) 



(5.9) 



On other other hand, applying the Burkholder-Davis-Gundy inequality, (2.11), and again using (2.13) we 
find that for < Ta < n < t A ^ A t'^P , 



E sup 

t'e[T^,Tt] 



JedW 



f-Tb \ 1/2 

<cE( / {d,a^{U),d,wfdt' 



<cE 



sup la.vi ( / {i + \\ur)dt 



<cE( /"Vv(t/)lL(u,y)|5.vpdt' 

1/2- 

< ^E 

e[ra,rt,] 



1/2 



(5.10) 



< ^E sup |a,v|2 +cE / "(1 + llC^f 
2 t'efT„,rhl 7t„, 



Gathering the estimates (5.7), (5.8), (5.9), (5.10) we apply the stochastic Gronwall Lemma, Proposi- 
tion 6.3 (from [GHZ09, Lemma 5.3]) and infer (5.3). This completes the proof of Proposition 5.2. □ 

5.2 L'^{M.) estimates for the Temperature and Salinity Equations 

We turn to the estimates for the temperature and the salinity, beginning with the estimates in L'^{M). 

Lemma 5.1. Suppose the standing conditions as in Theorem 2.1 are satisfied and let ([/, ^) he the maximal 
strong solution of (2.14). Define Tj^ as in (3.3). Then, for any t > 0, 



E sup (|r|i. + |5|i.) < 00, 



(5.11) 



Moreover, 



sup (|r|i4 + |S'|i4) < 00. a.s. 
t'e[o,tA{] 



(5.12) 



With defined as in (5.2), we have that VmiK^oa = 00. 
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Proof. We carry out the estimates for the temperature T. Those for the sahnity S are identicaL In order 
to find an evolution equation for |T|^4 we may proceed similarly to e.g. [Kry99, MROl]. We apply the 
Ito lemma to (2.1c) pointwise for almost every (x, z) G A4. This is justified given the regularity already 
established for T in Theorem 3.1. We then integrate the resulting system over M, apply the stochastic 
Fubini theorem (see [DPZ92]) and find that 



d\T\l,+ 12^T V / {d-jTfT^dM + 12vt f {d^TfT^^dM dt 

\ ,j JM JM j 

=4 / FTT^dMdt + &Y^ I aT{U)fT^dMdt + iY^ j (jT{U)iT^dMdWi 

JM />! ^ 1>1 



(5.13) 



Note also that we have once again used the cancelation properties of the nonlinear portion of (2.1c) (cf. 
(4.8)) in inferring (5.13). 

The external body forcing term is estimated as above (cf. (4.18)): 



FrT^dM 



M 



<c{\FT\% + l){\T\i, + l). 



(5.14) 



For the Ito correction term, using (2.13), wc find 



J2 I 'JT{U)fT'dM 

, JM 



<cJ2WTiU)i\UT\ 



<cWT{U)\l^^uy)\T\li 

<c(i + ||[/f)(i + |r|i.). 



(5.15) 



With an application of the BDG inequality, (2.11), wc find 

/ I] / (^T{U)iT^dMdWi <cK (^J^' J2 {^J '^T{U)iT^dAij dt 

<cEU^^'Y.\^T{u)i\UT\uit] 



E sup 



1/2 



<cE 



WmUyiy)\T\hdt 



1/2 



(5.16) 



<cE sup \T\lJ Til + \\Ur)\T\l,dt' 
<iE sup \T\l, + cE ril + \\Uf}{l + \T\Udt. 



Here, similarly to (5.10) above, Ta, Tb may be any stopping time such that < Ta < Tb < tATj^ A^. Combining 
the estimates (5.14), (5.15), (5.16) with (5.13) we apply the stochastic Gronwall lemma. Proposition 6.3, 
to infer (5.11). Now, due to Theorem 3.1, lim.K'fooT^ = oo. With this observation (5.12) and thus that 
limif-j-oo T^'^ ~ oo follow directly from (5.11). This completes the proof of Lemma 5.1 □ 

The final step is now to carry out the estimates for dzT and dzS in L^(A^): 
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Proposition 5.3. Assume the conditions of Theorem 2.1 consider (C/, ^) the maximal strong solution of 
(2.14). Take uk according to (5.2). Then, for any t > 0, 

E sup (|a,T|2 + |a,5n+ / {\dj:f + \\d,sf)dt'\<oo. (5.17) 

\_te[o,tAcrjcA4] Jo ) 

Moreover, up to a set of measure zero, 

sup (|a,T|2 + |9,5n+ / {\\d,Tf + \\d,Sf)dt' < 00. (5.18) 
te[o,tAC] Jo 

Proof. As above, in Lemma 5.1 we provide the estimates for T; tliose for S are identical. To tliis end we 
apply dz to (2.1c) and find: 

ddzT + a^[(v ■ V)r + w{w)dzT]dt - (nr'^^dzT + VTdzzdzT)dt = d^Frdt + a^crr(v, T, S)dW. (5.19) 

An application of the Ito formula then reveals: 

d|a,T|2 + 2\\dzTfdt =2(9,Ft, dzT)dt + |9,(tt(v, T, S)\l^(^^H)dt 

-2 / dz[{^ ■V)T + w{w)dzT]dzT dMdt + 2{dz(TT{^r,T,S),dzT)dW (5.20) 

JM 

■=iJi + JI + J3+ Jl)dt + JJdW 
For the nonlinear term, , an integration by parts reveals that, 

~Jl = y\ I {vjdjzT + dzVjdjT)d,TdM^y] [ {vjdjzTdzT ~ dzjVjTdzT ~ dzVjTd,jT)dM, 

J JM J JM 

and we may therefore estimate 

\Jl\ <c{\^\LA\d.T\\\dzT\L. + \\dM\\T\LAdzT\L. + l^.v^^ IT^^ ||a,T||) 

<c{\w\i^4dzT\\y^\dzT\^/^ + \\dM\\T\LAdzT\^/%dzTf/^^^^^ (5_2i) 

<c(|v|i.|9,Tp + iia.vf |T|i. + iia.vf + la.vpiTii,) + l(|a,rp + iia.Tf ). 

Also, 

Jj= / w{Y)dzTdzzTdM = -\ I dM^)idzT?dM^l f d^Vjid^Tf dM = - ( v^d.^Td^T dM, 

JM ^ JM ^ JM JM 



so that 



IM ^ J M ^ J M JM 

\JI\ < c\^,\L4^zn'^'\^zT\'^' < cHUdzTl' + ip.Tf . (5.22) 

The remaining drift terms are estimated directly using assumption (2.13): 

\Jl + Jl\ < c(l + \\Uf + \F\') + \\dzTf. (5.23) 

As above, in (5.10), we estimate the terms involving jj with the BDG inequality. We obtain for Ta < Tb < 
< A (Ta- A ^ 

E sup / Jr^dW 



< ^E sup IdzTl"^ + c£ f \l + \\U\\'^)dt' 



(5.24) 



With the stochastic Gronwall lemma. Proposition 6.3, the estimates (5.21), (5.22), (5.23), (5.24) imply 
(5.17). (5.18) now follows directly from (5.17), Proposition 4.1, Corollary 5.1 and Lemma 5.1. The proof of 
Proposition 5.3 is thus complete. □ 

Proof of Proposition 5.1 (conclusion). Due to Corollary 5.1 and Proposition 5.3 wc infer (5.4) and (5.18) 

(2) 

and hence (5.1). Now (5.1) implies that limR^oo = 00, completing the proof of the proposition. □ 
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6 Appendix: Auxilary Results 



We collect here, for the convenience of the reader, several technical results that have been used in an essential 
way in the analysis above. 

6.1 Pressure Estimates 

We have made use of a special case of [SvW86, Theorem 2.12] in the proof of Proposition 4.3. This result in 
[SvW86] provides L^U"^ estimates for the pressure terms appearing in the linear Stokes equation. 

Proposition 6.1. Suppose that d >2 and that A4o C M.'^ is a bounded open domain with smooth boundary. 
Assume that r G (1, 2) and that f e Lf^^{[0, oo); L'^{A4a)). If (q,p) solves the Stokes equation in M-q x [0, oo) 

dtq - J^Aq + Vp = f , V • q = 0, q\oMo = 0, 
then, for any < tq < ti < oo, 

r ivpii.(^„)ds < c ('iiq(To)ip + r \f\i.^Mo)A ' (6-1) 

J To \ J To J 

where c = c{d, r, A4o) is an absolute constant independent ofi, tq, ti. 

6.2 Two Nonstandard Gronwall Lemmas 

We made significant use of two non-standard versions of the Gronwall Lemma in the analysis above. The 
first version. Proposition 6.2, is used to close the estimates for v in for the proof of Proposition 4.1 
in Section 4. The second version. Proposition 6.3, has been useful for the analysis of stochastic evolution 
equation and is employed here in Theorem 3.2, Proposition 5.2, Lemma 5.1, and Proposition 5.3. 
The following proposition is somehow new: 

Proposition 6.2. Suppose that, for some t >0, we are given f,g,hG L^{[0,t]) and X G C([0,t)). Assume 
that X,f,g,h are all positive for a.e. t' G [0,t] and that f is continuous at t' = 0. //, for almost every 
< Ta < Tb < t, and for some fixed p > 1, we have 

sup X{t')<f{Taf+ sup X{t') f ' g{t')dt' + f \{t')dt', (6.2) 
then there exists an absolute constant c = c{t,p, |/(0)|, L^IlO such that, 

sup X(t') <c(l+ [ h(t')dt'\ . (6.3) 

t'G[0,t] \ Jo ) 

Proof. We begin with some preliminaries to determine the constant c in (6.3). Choose e > so that 
j^,g{t')dt' < 1/2 for any pair < s' < s < t such that s — s' < e. Now choose n in such a way that 
t/n < e/2. Then pick M > max{2, |/(0)|} large enough such that A(|/| > M) < t/^An), where A is the 
Lebesgue measure on the interval [0,t]. This later quantity M exists as a consequence of the Chebshev 
inequality and the fact that / G L^([0,<]). We now show (6.3) is satisfied by taking c = 2M^. 

To this end fix any t* G {t/2,t). Set to = 0, tn = t* . For each A; = l,...,n — Iwe may pick an element 
tk G {t*k/n,t*{k+ l)/n) so that f{tk) < M . Note also that according to this choice tk+i — tfc < e. Thus, for 
/s = 0, 1, . . . , n, by applying (6.2) with Tq = t^, n = tk+i, we infer, 

sup X{t') < il/P + i sup X(t') + f ^ h{t')dt'. 

t'e[tk:tk+l] ^ t'G[tk,tk+l] Jtk 
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Noting that X G C([0,i), we may rearrange and infer: 

sup X(t') < 2MP + 2 /" h{t')dt' < 2MP (l + [ h{t')dA . (6.4) 

The analogue of (6.3), with t* replacing t on the left hand side of the inequality now follows from a simple 
induction. Since M maybe chosen independently of t* G {t/2,t) the proof is now complete. □ 

The following 'stochastic Gronwall Lemma' has been established in [GHZ09]: 

Proposition 6.3. Fix t > and assume that X,Y,Z,R : [0,i) x O — > R are real valued, non-negative 
stochastic processes. Let t < t be a stopping time so that 



(6.5) 



E [ {RX + Z)ds < oo. 
Jq 

Assume, moreover that for some fixed constant k 

Rds < K, a.s. (6.6) 
Suppose that for all stopping times < Ta < < t 



J Y ds^ < coE (^X{Ta) + J [RX + Z) ds 



E sup X + I Yds] <co¥.[ X{Ta) + I {RX + Z)ds] , (6.7) 

where cq is a constant independent of the choice of Ta,Th. Then 

E\supX+[ Yds] <ce(x{0)+ [ Zds], (6.8) 
yte[o,T] Jq J \ Jo 

where c = c{cq, t, k). 
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